function [t, y] = f2(t, y, f)
f1 = f(t, y);
f2 = f(t + h/2, y + h*f1/2);
f3 = f(t + h/2, y + h*f2/2);
f4 = f(t + h, y+h*f3);
t = t+h;
y = h*(f1 + 2*f2 + 2*f3 + f4)/6;
    